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ABSTRACT 


The  motion  of  a  perfectly  conducting  fluid  past  a  fixed  body  in 
the  presence  of  an  aligned  magnetic  field  can  be  regarded  as  the  limit 
of  a  number  of  problems  viz:  (i)  when  the  fluid  is  a  finite  conductor, 
(ii)  an  unsteady  problem,  (iii)  the  magnetic  and  velocity  fields  are 
not  parallel  at  infinity.  The  limits  of  these  more  general  problems  agree 
in  predicting  a  force  on  the  body  and  on  the  existence  of  upstream  and 
downstream  wakes.  However,  if  the  magnetic  field  is  fairly  strong 
(1)  does  not  predict  a  downstream  wake  while  (il)  and  (iii)  do. 

An  attempt  to  reconcile  these  limit  results  is  made  here.  In  the 
first  part  we  show  that  the  unsteady  problem  associated  with  a  finitely 
conducting  fluid  and  an  aligned  field  is  non-unique  but  that  by  making 
an  appeal  to  the  theory  of  real  fluids  the  conclusion  of  (i)  and  (ii) 
can  be  obtained  as  special  cases  of  the  general  solution.  In  the  second 
part  the  steady  problem,  assuming  that  the  conductivity  is  finite  and  the 
magnetic  field  is  nearly  aligned,  is  considered  and  it  is  shown  that  the 
flow  fields  predicted  in  (i)  and  (iii)  are  limiting  cases  of  this  more 
general  problem.  The  manner  of  the  changeover  from  the  solution  for 
(i)  to  that  for  (iii)  is  elucidated. 


MAGNETO-FLUID  DYNAMICS  OF  BODIES  IN  ALIGNED  FIELDS 


K.  Stewartson 


1.  Introduction 

The  properties  of  the  two-dimensional  motion  of  an  incompressible 
inviscid  fluid,  of  electrical  conductivity  cr  ,  past  a  fixed  non-conducting 
body  end  in  the  presence  of  an  imposed  magnetic  field  have  attracted  much 
interest  in  recent  years.  In  particular,  the  properties  of  the  motion  when 
the  magnetic  and  velocity  fields  are  aligned,  i.  e.  parallel,  at  large  distances 
from  the  body  and  o-  is  infinite  have  been  the  subject  of  some  controversy. 
Sears  and  Resier  (1)  pointed  out  that  in  these  circumstances  there  is  a 
simple  solution  of  the  equations  of  motion  in  which  the  magnetic  field 
outside  the  body  is  everywhere  parallel  to  the  velocity  field  which  in  turn 
is  identical  with  that  when  o-  =  0 .  Inside  the  body  the  magnetic  field 
vanishes.  One  Important  deduction  from  this  solution  is  that  the  drag  on 
the  body  is  zero.  They  also  briefly  discussed  the  solution  when  o-  is 
large  but  not  finite  and  concluded  that  if  the  body  is  thin,  their  solution 
is  consistent  with  it  in  a  limited  sense.  Subsequently,  Ring  ( 2)  extended 
their  solution  to  include  the  effects  of  oscillating  a  thin  body  about  an  axis. 
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Some  doubts  about  the  relevance  of  these  solutions  to  physically 

realisable  situations  were  cast  by  Stewartson  (3)  ,  who  meant  that  It  was 

doubtful  whether  the  solution  obtained  by  Sears  and  Rosier  (1)  was  the 

limit  of  a  general  fluid  dynamic  problem  in  which  the  motion  was  started 

started  from  rest,  at  which  time  the  magnetic  field  was  everywhere  uniform 

and  at  an  inclination  a  to  the  mean  direction  of  flow,  and  the  fluid  was  of 

finite  conductivity,  as  t  -►  <»,  a  -•  0 ,  <r  •♦  <»,  either  jointly  or  separately. 

In  support  of  these  doubts  he  considered  the  unsteady  motion  of  a  thin  body 

with  (T  =  <»,  a  =  0 ,  showing  that  the  disturbance  produced  manifested  itself 

partly  as  a  potential  flow  and  partly  In  the  form  of  wakes  travelling  from  the 

body  along  the  lines  of  magnetic  force  without  being  either  dispersed  or 
2  2 

diffused.  Let  P=H_/4irpV^,  where  H^,  are  the  magnitudes  of  the 
undisturbed  magnetic  and  velocity  fields  and  p  is  the  fluid  density:  then 
if  p  >  1  these  wakes  travel  upstream  and  downstream  of  the  body  while  if 
P  <  1  both  travel  downstream-  Forthe  simple  solution  tobecorrect  these 
wakes  must  be  absent  but  all  that  could  be  deduced  with  certainty  from  the 
temporal  solution  is  that  the  assumption  of  small  disturbances  breaks  down 
at  infinite  times.  It  Is  noted  however  that  the  argument  in  the  paper  is  not 
entirely  satisfactory.  At  any  finite  value  of  t  It  may  be  shown  that  the 
magnetic  field  is  continuous  at  the  body  ( taking  the  permeability  of  both  fluid 
and  body  to  be  unity)  because  the  normal  component  of  the  magnetic  field 
is  non-zero.  In  the  paper  it  was  tacitly  assumed  that  the  same  condition 
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must  hold  at  an  Infinite  time  but  as  M.  B.  Glauert  kindly  pointed  out  to 
the  author  the  properties  of  the  steady  state  solution  imply  that  the  normal 
component  of  the  magnetic  field  must  be  zero  if  the  disturbances  are  small 
so  that  a  discontinuity  in  the  tangential  component  of  the  magnetic  field  is 
permissive.  If  the  discontinuity  is  allowed  and  the  disturisances  in  the  fluid 
are  kept  small  the  solution  obtained  by  Sears  and  Resler  ( 1)  is  recovered, 
but  of  course  the  disturbance  to  the  magnetic  field  in  the  body  is  not  small 
since  it  is  actually  zero.  Consequently,  as  stated  above,  the  only  legitimate 
conclusion  from  a  study  of  the  general  temporal  problem  is  that  the  ultimate 
steady  solution,  if  it  exists,  is  not  uniformly  a  small  perturbation  from  the 
undisturbed  state. 

There  is  however  a  limiting  temporal  solution  which  can  be  carried 
out  without  any  such  difficulty,  namely  if  p  »  1  and  here  Stewartson  (3) 
showed  that  the  ultimate  solution  is  quite  different  from  the  prediction  in 
(1)  and  contains  two  wakes  extending  upstream  and  downstream  to  infinity. 
Thus  at  least  one  of  these  two  solutions  must  be  a  singular  limit.  If  infinite 
wakes  are  permitted  it  follows  that  arbitrary  conditions  at  infinity  cannot  be 
Imposed  on  the  flow  when  p  >  1  without  losing  physical  reality  and  the 
question  of  determining  physically  realistic  solutions  becomes  an  order  of 
magnitude  harder  to  answer.  We  note  that  the  solution  with  p  »  1  is 
obtained  by  neglecting  the  convection  terms  in  the  equations,  1.  e.  by 
assuming  that  the  fluid  has  velocity  but  does  not  move! 
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The  second  limit  problem  a  -»  0  was  investigated  by  Stewartson 
(4)  assuming  a  prefectly  conducting  tluio,  and  steauy  conoitions.  In 
fact  the  general  problem  when  the  magnetic  field  is  at  an  arbitrary  angle  of 
incidence  a  to  the  velocity  field  was  solved  and  on  taking  the  limit  or  0 
it  was  found  that  the  flow  properties  axe  in  accord  with  those  suggested  by 
the  temporal  solution  and  not  with  the  hypothesis  that  wakes  are  absent. 

In  particular  the  drag  on  a  symmetric  and  symmetrically  disposed  body 
takes  the  form 

4pV^  y 

— ^/{f‘(x)}  dx  (1.1) 

where  f''(x)  is  the  slope  of  one  side  of  the  body;  if  the  flow  properties 
were  Independent  of  the  magnetic  field  the  drag  would  be  zero.  However, 
this  argument  is  not  conclusive  since  the  theory  breaks  down  in  the  limit 
a  -»  0  because  the  disturbances  are  no  longer  small  and  because  when 
a  <<  (  ,  where  c  is  a  measure  of  the  slope  of  the  body,  the  wakes 
Intersect  the  body. 

The  third  limit  problem  is  to  study  the  steady  flow  properties 

when  a  =  0  but  v  is  finite  with  particular  reference  to  the  solution 

when  <r  is  large.  This  has  been  carried  out  by  Lary  (5)  for  a  thin 

1 

body  who  found  that  the  drag  is  proportional  to  ir  ‘  when  r  is  large. 

Further  if  P  <  1  a  strong  wake  develops  downstream  of  the  body  only, 
as  (T  is  increased,  while  if  P  >  1  the  strong  wake  develops  upstream  only. 

Since  the  disturbances  produced  by  the  body  Increases  with  cr  the  solution 
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obtained  by  Lary  (5)  is  of  limited  validity,  the  precise  condition  being  that 
2 

(T  c  <  <  1  ,  and  so  cannot  be  applied  immediately  to  the  case  or  =  <» . 
However  it  does  present  further  evidence  in  support  of  the  view  th  at  the 
flow  of  a  fluid  of  infinite  conductivity  is  grossly  different  from  that  of  a 
fluid  of  zero  conductivity,  that  it  contains  either  downstream  or  upstream 
wakes  and  that  the  force  on  the  body  is  non- zero.  Nevertheless  the 
solutions  obtained  by  the  three  limiting  procedures  t-*«o,  a-*0,  (r-*oo 
do  differ  in  certain  details  of  which  the  most  important  is  that  in  Lary's 
solution  there  is  no  downstream  wake  when  P  >  1  while  in  the  other 
limiting  procedures  (t  -►  <»  and  a  -*  0)  such  a  weike  is  present.  In  this 
paper  we  shall  attempt  to  reconcile  the  limiting  procedure  o-  -►  »  with  the 
other  two.  First  we  shall  discuss  the  properties  of  the  unsteady  linearized 
equations  in  the  presence  of  an  aligned  magnetic  field  and  show  that  the 
equations  do  not  have  a  unique  solution  unless  appeal  is  made  to  the  theory 
of  real  fluids  in  order  to  determine  the  vorticity  at  the  rear  of  the  body.  It 
is  argued  that  if  the  flow  remains  almost  uniform  and  the  body  is  appropriately 
streamlined  the  non-uniqueness  may  be  removed  and  Lary's  solution  is  then 
recovered.  In  the  temporal  problem  with  v  »  1  for  which  a  solution  could 
be  found  p  »  1  the  streamlines  are  greatly  disturbed  and  there  are  some 
theoretical  grounds  for  accepting  the  solution  in  (3)  if  o-  =  and  in  (4) 


if  (T  <  00  . 


-6- 


#360 


In  the  second  problem  the  steady  flow  properties  when  or  »  1  , 

2 

a  «  1  but  0-  a  is  finite  are  studied  and  the  manner  of  the  transformation 

of  the  solution  from  (5)  to  (6)  as  varies  from  to  0  is 

demonstrated.  In  particular  the  disappearance  of  the  strong  wake  down¬ 
stream  of  the  body  as  <rcr^  -»  0  for  §  >  1  is  elucidated. 

Although  it  seems  clear  that  the  simple  solution  obtained  by  Sears 
and  Resler  (1)  is  not  the  limit  of  any  of  the  more  general  problems  considered 
here,  we  are  not  yet  justified  in  concluding  that  it  is  essentially  isolated.  For 
there  is  another  limit  which  has  not  yet  been  considered,  namely  v  0  ,  where 

V  is  the  kinematic  viscosity  of  the  fluid.  In  the  discussion  it  has  been  tacitly 
assumed  that  v  =  0  from  the  outset  and  it  may  be  that  if  the  order  of  taking 
limits  is  reversed  the  flows  obtained  may  be  quite  different.  In  this  connection 
it  is  noted  that  Hasimoto  (6)  has  shown  that  if  v=<»,  v^O,  |3<1,  there 
is  a  solution  of  the  governing  equations  which  is  identical  with  that  for  a  non¬ 
conducting  fluid  of  different  viscosity:  if  p  >  1  the  two  flows  are  closely 
related.  Whether  Hasimoto's  solution  is  the  limit  of  solutions  with  finite  o-  , 
as  (r-»  00,  is  an  unresolved  question  at  present  but  certainly  its  limit  as 

V  -•  0  is  identical  with  (1)  . 
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PAKT  I 

UNSTEADY  THEORY 


2,  The  Basic  Equations 

Consider  a  cylinder  body  of  characteristic  length  a  fixed  in 
an  incompressible  fluid  of  electrical  conductivity  v  and  permeability 
t*  =  1.  The  fluid  is  initially  at  rest  and  at  time  t'  =  0  is  set  in  motion 
such  that  the  velocity  at  an  infinite  distance  upstream  of  the  body  changes 
abruptly  to  a  magnitude  in  a  direction  perpendicular  to  the  generators 
of  the  body  (i.  e.  to  the  z  direction)  and  subsequently  is  kept  constant. 

In  addition  a  magnetic  field  is  imposed  which  at  all  finite  times  is  uniform, 
of  magnitude  and  parallel  to  the  velocity  field  at  infinite  distances 
upstream  of  the  body.  Denoting  the  velocity  and  magnetic  vectors  by 
the  fluid  pressure  by  p  and  the  density  by  p  the  governing 
equations  are  . 

1 

^  -  curl  ^  curl  curl  H  ^  0  ,  dlv  H  =  0  ( 2. 1 ) 


at' 


+  (a;grad)a,= -- grad  p +  — (curl  IIxH),  div  a,  =  0  .  (2.2) 


Let  ax,  a^  measure  distance  along  and  perpendicular  to  the  direction  of  the 
undisturbed  stream  and  in  the  plane  of  motion;  let  at'/V^  =  t .  Further  suppose 
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that  the  disturbances  caused  by  the  body  are  small  so  that  squares  and 
products  of  the  perturbations  in  q  ,  H  may  be  neglected.  Write 

^  =  =  +  (2.3) 

where  ^  is  a  unit  vector  In  the  x  direction  and  on  substituting  into  ( 2. 1), 

( 2.  2)  we  get 


^  -  R 


at”  8x“  “  ®  *  8t”  ■*'  ^ 


8Jl 

gradP  +  p—  (2.4) 


where 


R  =  4ira«rV^,  P  =  4iTpV^  , 


(2.5) 


p  =  -pV^(P  +  Ph^)  +  p^ 


(2.6) 


and  p^  is  the  fluid  pressure  at  infinity.  The  velocity  and  magnetic  fields 
may  be  expressed  in  terms  of  two  scalar  functions  ^  as  follows: 


pl!lL 

8x8y 


_  8  »  8  » 

''y  ■  8x8y  8y8t 


8x 


(2.7) 


sfl  8^ 

X  „  2  ^  8y8x  8y8t 
ox 


2  2  2 

y  "  8xay  ’  -  2  "  8x8t 
ox 


(2.8) 
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where 


and 


®!i+sjfe  =0 

8x2  3^2 


(2.9) 


( 2. 10) 


Equation  (2.10)  Is  standard  but  (2.11)  Is  unusual  and  has  some  interesting 
features.  If  t)<  is  Independent  of  ^  it  reduces  to 

—  [  R(p.l)|i]  =0  (2.12) 

of  which  one  solution  is  independent  of  x  and  the  others  satisfy  an  equation 
reminiscent  of  the  vortlclty  equation  in  Ose6n's  approximation  to  the  Navier 
Stokes  equations  of  classical  hydrodynamics.  The  existence  of  a  solution 
independent  of  x  means  that  there  is  a  possibility  that  a  wake  may  develop 
ultimately  extending  to  infinity ,  which  in  turn  may  mean  that  it  is  not  uniformly 
valid  to  Impose  uniform  conditions  at  infinity.  Let  us  examine  how  this  solution 
arises  by  looking  at  the  variation  of  with  t .  It  is  convenient  to  take  the 
Laplace  transform  of  (2.11)  with  respect  to  t,  using  £  as  new  parameter 
and  denoting  the  transformed  function  by  4*  •  The  equation  satisfied  by  ^  is 
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+  =0  (2.13) 

of  which  solutions  can  be  found  by  writing 

+  =  /"eik.  dk  (2.U) 

0  ^ 

where  k  is  real  and  6  such  that  the  real  part  of  Ox  is  negative.  From 
(2.13)  6  satisfies 

(0  +  s)(e^-  kS  -  R(e  +  s)^+  Rpe^  =  o  (2.15) 

and  we  take  £  to  be  real  and  positive.  When  6  =  -eo  the  left  hand  side  of 

2 

(2.15)  is  negative,  when  6  = -s  it  is  positive  and  equal  to  RPs  ,  when 

2  2 

0=0  it  is  negative  and  equal  to  -(sk  Rs  )  and  finally  when  0  =  +<»  it 
is  positive.  Thus  there  are  three  real  roots  0^,  0^,  0^  of  (2.15)  which 
satisfy 

0j  <  -s  <  02  <  0  <  0j  (2.16) 

the  equality  signs  being  possible  only  if  8=0,  P  =  or  R  =  »,  or  if 
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p  s  0  and  8  =  k .  Otherwise  the  roots  are  distinct.  Further  of  the  roots 
only  0^  Is  relevant  to  large  positive  values  of  x  and  8^ ,  0^  to  large 
negative  values  of  x ,  since  «)«  cannot  be  exfxjnentlally  large  at  Infinity. 

When  s  Is  large  ( 2. 16)  reduces  to 

(8+  s)(0^-  Rs)  =  0  (2.17) 

so  that  0j  cs  -s  ,  ®2  “  ■  »  0^  —  +  as  s  -►  «> .  These  roots 

may  be  Interpreted  In  terms  of  solutions  of  the  equation  for  .  If  0  =  -s 
when  s  Is  large  It  follows  from  ( 2. 14)  that 

4<aiF(x-t,  y)  (2.18) 

when  X  - 1  Is  small,  where  F  Is  some  function  at  present  arbitrary,  and 
corresponds  to  a  disturbance  travelling  downstream  with  the  main  stream 
velocity.  Such  a  solution  might  have  been  anticipated  since  the  linearization 
of  the  equation  Is  partly  equivalent  to  assuming  that  convection  Is  taking 
place  with  the  main  stream  velocity  rather  than  with  the  local  stream  velocity. 

If  0  =  -  ''TSa  It  follows  from  ( 2. 14)  that 


'|/~g(x,  y,  t)e 


-Rx^/4t 


(2. 19) 
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when  t  is  small,  where  g  is  also  arbitrary.  This  solution  corresponds  to 
a  diffusion  and  arises  if  on  setting  the  fluid  into  motion  the  change  in  the 
magnetic  field  due  to  the  currents  engendered  leads  to  an  initial  discontinuity 
in  either  of  its  components  at  the  surface  of  the  body.  The  discontinuity  begins 
to  diffuse  in  a  manner  characterized  by  ( 2. 19) . 

We  now  enquire  what  happens  to  these  solutions  as  t  «  assuming 
that  the  ultimate  motion  is  steady.  The  assumption  is  equivalent  to  assuming 
that  6i,  B  are  regular  functions  of  £  in  the  half  plane  re^>0  and  that  on 
the  imaginary  axis  of  s  the  singularities, if  any, are  sufficiently  weak,  except 
at  s  =  0 ,  for  their  contribution  to  t)^  to  die  out  algebraically  as  t  .  It 
is  then  sufficient  to  consider  the  behaviour  of  6  near  s  s  0.  When  s  =  0 
( 2. 1  )  reduces  to 


0(0^-  k^)+  R(p-1)0^  =  0  ; 


( 2.  20) 


from  the  definition  of  0^,  6^,  0^  in  (2.16)  it  follows  that 


0j  --^R(p-l)  -^^/R^P-l)^  +  4k^ 


02-0 


( 2.  21) 


03->7R(P-1)  +  ^'>/R^P-1)^  + 4k^ 


J 
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as  s  **  0 .  More  precisely  as  s  0 

e^/s--!  (2.22) 

One  may  reasonably  suppose  that  the  correct  form  for  vt<  in  a  particular 
problem  Is  compounded  of  multiples  of  the  three  solutions  defined  by 
( 2. 14)  ( 2. 15)  and  further  that  that  part  of  the  disturbance  characterized 
by  v|i  is  centred  at  the  body  at  t  =  0 ,  moving  away  from  It  as  t  Increases. 

It  then  follows  that  only  the  solution  i|i^  given  by  (2. 13)  with  6  =  0^  affects 
the  region  upstream  of  the  body  (where  x  is  negative),  that  when  t  is  small 
Is  diffusive  and  that  when  t  is  large  4*^  decays  exponentially  with  x  for 
finite  R  and  p  #  1  .  Hence  at  large  distances  upstream  the  perturt>ation  in 
the  state  of  the  fluid  is  dominated  by  the  harmonic  term  ^  .  Downstream  of 
the  body  the  solutions  with  6  =  6^,  6^  are  relevant  and  a  comparison 

between  ( 2. 17)  ( 2.  21)  ( 2.  22)  shows  their  roles  interchange  as  t  increases 
from  0  to  When  t  is  small  4^^  is  wave -like  while  is  diffusive. 

On  the  other  hand  when  t  is  large  decays  exponentially  with  x  while 
4*2  is  s  function  of  x-t  and  y.  In  fact  since  we  are  assuming  steady 
motion  we  have  from  (2.4),  (2.8)  that  is  independent  of  x 

and  4*2  has  the  form 

(x-t)F(y). 

Of  the  solutions  derived  4/^ ,  41^  also  satisfy 
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V\  +  R(P-1)|J  =0  (2.23) 

and  for  convenience  we  shall  denote  the  x-derivative  of  that  part  of  t|/  which 
satisfies  (2.23)  by  'i'  .  The  other  solution  from  continuity,  only 

zero  inside  a  region  C  bounded  by  straight  lines  parallel  to  the  x  axis 
and  by  Infinity  on  the  downstream  side.  The  fourth  part  of  the  boundary  can 
be  defined  when  the  body  is  thin:  it  conslstc.  of  the  rear  portion  of  the  body 
and  the  straight  lines  mentioned  above  touch  the  body. 

If  the  body  is  not  thin  the  linearized  equations  are  in  general  not  valid 
in  the  neighborhood  of  the  body,  so  that  C  cannot  be  fully  defined.  However 
the  description  of  \\t  as  t  given  above  can  be  used  to  discuss  the  ultimate 
flow  pattern  at  large  distances  from  the  body.  It  will  consist  of  a  harmonic 
term  and  a  downstream  wake  for  all  values  of  R,  p  .  In  addition  if  R  is 
large  there  is  a  wake -like  contribution  from  ^  extending  downstream  in  a 
region  somewhat  like  C  if  p  <  1  and  extending  upstream  in  an  otherwise 
similar  region  if  p  >  1  .  In  fact  if  R  =  <»,  P  <  1 ,  4'  is  of  the  same  form  as 
Sij^^/dx  when  x  is  large  and  positive  and  if  p  >  1,  'Jr  is  Independent  of  x 
when  X  is  large  and  negative.  It  is  noted  that  makes  no  contribution 
to  ^  in  virtue  of  (  2.  8)  so  that  at  finite  R,  p  only  the  velocity  distribution 
can  be  non-uiUform  in  C .  On  the  other  hand  ^  can  be  non-zero  in  the  wakes 
produced  by  '1'  when  R  =  . 

Let  us  now  examine  the  properties  of  tjj  in  the  two  limiting  cases 
R-*  0  and  R-*<».  First  if  R  =  0  (2.15)  reduces  to 
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(6+  s)(e^-  k^)  =  0  (2.  24) 

with  roots  U  =  -s  and  6  :fck  and  here  the  root  which  behaves  like  6^ 
when  s  Is  large  behaves  like  6^  when  s  is  small  and  vice-versa.  There 
Is  no  interchange  of  roots  at  s  =  k  when  ( 2.  24)  has  a  double  root  as  may 
be  seen  by  tracing  the  roots  when  arg  s  is  small  and  positive.  Thus  the 
solution  in  the  limit  R  -*  0  4  is  not  the  same  as  the  solution  when  R  =  0 
and  so  care  must  be  exercised  in  inferring  flow  properties  when  R  is  small 
from  those  when  R  =  0  .  Sketches  of  the  behaviour  of  6^ ,  6^  regarded  as 
functions  of  s  when  R  =  0  and  when  0  <  R  «  1  are  given  in  Figure  1  . 


FIGURE  1 
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Second  if  R  =  <»  (2.15)  reduces  to 
(0  +  s)^  =pe^ 

1.  e. 

6  =  -1)"^  . 

These  two  solutions  correspond  to  6^  6^  if  P  <  1  and  to  6^  >  0^  if 
p  >  1  .  The  third  root  of  (2. 15)  becomes  +«  if  P  <  1  and  -»  if  p  >  1  . 

On  the  other  hand  if  R  is  large  but  finite  *♦  -1  as  s  -•  0  which  is 
different  from  ( 2.  25) .  A  more  careful  scrutiny  of  ( 2. 15)  shows  that  when 
R  is  large  and  s  is  small  the  form  taken  by  the  roots  depends  crucially  on 
Rs.  Take  P<1  and  then  0j=R(l-p)  which  is  effectively  infinite.  The 
other  two  roots  satisfy 

-k^(0  +  s)  =  R(e  +  s)^  -  RP0^  (2.26) 

with  solution 


(2.25) 


0  k^-t-  2R8*  N/k^  4pk^Rs  4P(Rs)^ 

s  "  2(l-p)Rs 


(2.  27) 


Using  the  inversion  formula  for  Laplace  transforms  we  can  now 
describe  the  behaviour  of  ^  when  t  is  large.  Ifp<l,  l«t«R, 
the  solutions  of  (2.13)  corresponding  to  0^,  0^  are  two  Alfv6n  waves 
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travelllng  downstream  from  the  body  without  dispersion  and  with  velocities 
given  by  1  ^  .  As  t  increases  through  R  dispersion  and  change  of 

speed  of  propagation  occur  until  when  t  »R  the  flow  pattern  described 
in  the  general  case  is  recoveibd.  Again  therefore  caution  must  be  exercised 
in  inferring  steady  flow  properties  when  R  is  large  from  those  when  R  is 
infinite.  This  point  has  already  been  discussed  by  the  author  in  a  particular 
case  (7)  to  which  the  reader  is  referred  for  further  details. 
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3.  Application  To  Flow  Past  Thin  Bodies 

The  linearized  theory  of  magneto-fluid  dynamics  may  be  used  to  study 
the  complete  flow  field  in  two  circumstances:  when  the  body  is  thin  and 
orientated  so  that  the  streamlines  are  almost  straight  and  when  the  magnetic 
field  is  strong  (P  »  1)  so  that  the  Inertia  terms  are  negligible.  In  this 
section  we  consider  the  first  of  these  deferring  the  second  to  §4  .  Suppose 
the  body  occupies  the  neighborhood  of  that  part  of  the  x  axis  for  which 
|x|<l.  From  the  theory  of  §2  it  follows,  on  assuming  that  the  ultimate 
flow  field  is  steady  and  almost  uniform,  that  in  the  fluid 


V 

X 


8x  8y 


+ 


F(y)  , 


V 

y  8y  8x 


h 


X 


a»  8* 
8x  8y 


h  -  £?  .  13! 

y  “  8y  ”  8x 


(31) 


where  we  have  written  for  8^/8x  ,  »  satisfies  (2.23)  and  both 
grad  4,  grad  ^  vanish  at  Inflnityt  the  box  surrounding  F(y)  in  the 
equation  for  v  means  that  there  is  an  additional  contribution  to  v  , 
independent  of  x ,  inside  ,  a  region  bounded  by  the  rear  part  of 
the  body  and  two  lines  parallel  to  the  x  axis  which  touch  the  body  and 
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extend  to  infinity  downstream  of  it  ( See  Figure  Z) . 


FIGURE  2 


Inside  the  body  the  velocity  field  is  zero  and  the  velocity  field  is 

harmonic.  As  a  consequence  the  variation  in  the  magnetic  field  across  the 

body  is  small  and  since  the  tangential  component  of  the  magnetic  field  and 

the  normal  component  of  the  Induction  are  continuous  at  the  surface  of  the 

body,  it  follows  that  for  fixed  x  the  values  of  h  ,  h  on  the  top  and 

X  y 

bottom  surfaces  of  the  body  are  equal.  A  further  condition  to  be  satisfied 
on  the  body  is  that  the  velocity  be  tangential  to  the  surface:  if  the  surfaces 
of  the  body  are  given  by  y  =  f^{x)  ,  the  *  signs  denoting  the  top  and  bottom, 
this  condition  reduces  to 


V 

y 


=  f;  (X) 


(3.  2) 


on  the  body,  effectively  y  =  0± ,  |x|  <  1  .  There  is  now  sufficient 
Information  formally  to  complete  the  determination  of  *  ,  ijr  and  this 
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has  been  done  by  Lary  (5)  In  a  sufficiently  representative  number  of  cases. 

Lary  tacitly  assumed  that  F(y)  ^  0  but  so  long  as  it  is  small  it  can  be 

regarded  as  a  quantity  to  be  found  after  ^ ^  are  known  and  in  particular 

makes  no  contribution  to  the  forces  on  the  body.  Hence  his  description  of 

these  forces  Is  essentially  complete.  For  full  details  of  his  solution  the 

reader  is  referred  to  his  paper  but  it  is  noted  that  for  large  R  a  strong 

wake  develops  upstream  or  downstream  according  as  P  >  1  or  P  <  1  and 
1 

that  the  drag  oc  R  >  .  This  means  that  the  theory  becomes  invalid  for 

sufficiently  large  R  and  he  showed  that  the  condition  for  validity  is  that 
-2 

R  «  c  where  c  is  a  representative  thickness  of  the  body.  There  remains 
the  determination  of  F(y)  which  since  it  leads  to  a  downstream  wake  for  all 
p  is  the  essential  link  between  the  solutions  with  c  «  1  and  those  with 
p  »  1  . 

As  we  have  seen, the  presence  of  F(y)  is  bound  up  with  the  root 
6  =  62  of  the  cubic  equation  ( 2. 14)  and  corresponds  to  magnetic  diffusion 
from  the  surface  of  the  body  at  t  =  0 .  Thus  the  solution  corresponding  to 
6  =  62  is  non-zero  at  any  finite  time  and  a  further  argument  is  necessary 
at  t  =  00  before  it  can  be  included.  Basing  the  argument  on  uniform 
conditions  at  infinity  is  unsatisfactory  since  the  effect  upstream  of 

e 

improving  such  a  condition  at  large  positive  values  of  x  dies  out 
algebraically  in  ^  and  exponentially  in  ^  so  that  it  cannot  affect  the 
flow  near  the  body. 


#360 


-21- 


The  further  condition  by  which  F(y)  is  determined  must  be  sought 
in  the  flow  properties  near  the  body  and  we  suggest  that  the  clue  to  the 
condition  lies  in  the  vortlclty.  The  necessity  for  a  condition  on  the  vortlclty 
in  a  linearized  theory  is  most  easily  seen  from  a  study  of  the  flow  when  the 
conductivity  of  the  fluid  is  zero.  The  governing  equations  are  then 

^  ^  =  grad  P,  dlv  ^  =  0  (4.  3) 

from  (  2.  5) ,  (  2.  6) .  In  two-dimensional  motion  (4.  3)  implies  that 


at  8x 


(4.4) 


where  ^  is  the  vortlclty.  This  equation  simply  means  that  vortlclty  is 
convected  along  lines  parallel  to  the  x  -  axis .  On  as  suming  that  all  vortlclty 
originates  at  the  body  it  follows  from  (4.  4)  that  ^  is  zero  except  inside 
,  which  extends  downstream  from  the  body.  Inside  the  behaviour 
of  ^  cannot  be  found  without  making  some  assumption  about  ^  on  that 
part  of  the  body  bounding  {ABC  in  Figure  2  }  .  Thus  if  ^  ss  g{y,  t)  there, 


4  =  8{y,  t  -  X  +  X(y)) 

in  where  x=X(y)  on  ABC;  as  t » 


4-*g(y,  00) 


(4.5) 
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for  all  X .  When  the  conductivity  of  the  fluid  is  zero  the  additional  assumption, 
by  which  g  may  be  found,  emerges  from  a  study  of  the  flow  properties  of  a 
fluid  with  small  viscosity  as  follows.  It  is  supposed  that  the  shape  of  the 
body  is  such  that  the  boundary  layer  rei..ains  attached  to  it  over  almost  all 
of  its  length.  Then  over  almost  the  ^ole  width  of  the  fluid  particles 
have  zero  vorticity  by  Kelvin's  circulation  theorem  since  they  had  none 
upstream.  Hence  g  =  0  except  possibly  when  y/Ac  «  1 . 

If  the  fluid  is  conducting  ( 4. 4)  becomes 


dt 


a,  a  8h  6h 

ii  =  6-5-  SI  - 

8x  '^Sx  8x  8y 


(4.6) 


and  again  i  can  be  given  an  arbitrary  value  on  ABC .  In  order  to  fix  t,  and 

obtain  a  well-posed  problem  we  adopt  a  similar  argument  to  that  for  o-  =  0  . 

First  we  suppose  that  the  body  is  streamlined  so  that  an  Insignificant  amount 

of  vorticity  is  shed  from  the  boundary  layer.  Second  we  make  use  of  the 

magneto-fluid  dynamic  extension  of  Kelvin's  theorem  in  a  linearized  form. 

This  implies  that  the  vorticity  at  may  be  calculated  from  the  vorticity 

at  whose  y-coordlnate  is  the  same  as  that  of  (see  Figure  2)  by 

integrating  (4. 6)  from  Dj  to  using  the  values  of  h^,  h^  on  the 

boundary  of  the  body,  i.  e.  by  the  values  of  h  ,  h  as  the  line  y  =  0 . 

X  y 

Here  the  motion  is  steady  and  so  the  difference  in  the  values  of  the 
vorticity  at  and  is  equal  to 
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Hence  on  comparing  the  values  of  the  voiticity  at  , 

F{y)  »  0 

and  Lary's  solution  is  completely  valid  if  R«  c"^  . 
passing  that  if  the  vorticity  condition  had  led  to  a  non 
it  would  have  been  0(€^)  since  F*(y)  =0(€) 


it  follows  that 

It  is  worth  noting  in 
-zero  value  of  F(y) 
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4,  Discussion 

It  was  pointed  out  in  §  1  that  the  main  diffeieuce  between  the  unsteady 
and  the  finite  conductivity  approaches  to  the  flow  of  a  perfectly  conducting 
fluid  past  a  body  in  the  presence  of  an  aligned  field  occurs  if  p  >  1  and 
concerns  the  existence  of  a  velocity  wake  on  the  downstream  side  of  the 
body.  It  has  been  shown  in  §3  that  the  steady  state  solution  of  the  problem 
when  0-  <  0°  concalns  an  arbitrary  function  which  gives  rise  to  this  wake  and 
is  associated  with  the  vorticlty  shed  frxxn  the  body.  Provided  the  flow  remains 
almost  uniform  and  the  body  is  streamlined  to  delay  boundary  layer  separation 
the  vorticlty  shed  may  be  neglected  and  this  function  vanishes.  The  flow 
pattern  is  then  as  predicted  by  Lary.  Thus  the  uncertainty  in  the  solution 
is  resolved^  by  an  appeal  to  the  realism  of  the  flow  predicted. 

Unfortunately  the  flow  diverges  strongly  from  its  uniform  state  either 
asR-»<»  or  IfRsw  and  p  finite  there  is  a  controversy  about 

the  character  of  the  flow  but  for  p  =  <»,  R  =  <»  there  exist  strong  magnetic 
and  velocity  wakes  upstream  and  downstream  of  the  body  {i).  Although  the 
solution  with  R  finite  and  p  =  «>  (7)  differs  from  that  one  because  the 
magnetic  wake  is  absent  the  two  can  be  reconciled.  Consequently  we  need 
only  consider  the  discrepancy  between  the  solutions  with  p  »  1  and  finite 
p  ,  R  being  finite  in  each  case. 

Let  us  now  examine  the  form  taken  by  ( 3. 1)  when  p  is  large  and  R 
finite.  The  equation  satisfied  by  reduces  to 


#360 


-25- 


8x 


=  0 


(4.1) 


except  near  the  part  ABC  of  the  body  where  a  boundary  layer  type  of 
phenomenon  can  develop  In  St .  Assuming  that  the  velocity  Is  of  order 
one  at  most,  h^ ,  h^ ,  ,  •  are  negligible  In  ( 4. 1)  and 


Fj(y) 

In 

F^Cy) 

In  C_ 

(4.  2) 

0  elsewhere 

where  F^,  F2  are  functions  of  y  to  be  found.  The  solution  for  P»1 
(7)  Is  obtained  by  setting  F^  s  F2  =  -1  so  that  the  velocity  Is  zero  In 
and  undisturbed  everywhere  else.  This  form  Is  consistent  with  the 
approximations  made  and  appears  to  be  more  realistic  than  the  form 
obtained  by  extending  Lary's  condition  on  F(y)  for  then.  »  0  and 
F2  =  -1  :  the  boundary  condition  on  the  velocity  over  ABC  being  satisfied 
by  means  of  the  boundary  layer  In  St  mentioned  earlier.  The  alternative 
requires  motion  of  the  fluid  In  ABC  which  a  real  fluid  would  be  unlikely 
to  sustain. 

There  is  however  some  experimental  evidence  which  should  be 
borne  In  mind  in  connection  with  these  problems.  The  theoretical  flow 
properties  when  p  »  1  are  very  similar  to  those  when  a  body  moves 
slowly  along  the  axis  of  a  rotating  fluid  namely  that  upstream  and  down¬ 
stream  wakes  are  predicted.  Experiments  to  test  this  theory  have  been 
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carried  out  by  Taylor  (8)  who  examined  the  flow  in  the  upstream  wake  and 
by  Long  ( 9)  who  examined  the  flow  ahead  and  behind  the  body.  Both 
experiments  confirmed  the  existence  of  the  upstream  wake  but  Long  found 
that  the  fluid  streamlines  closed  up  behind  the  body  contrary  to  the  theoretical 
prediction.  No  firm  explanation  of  the  discrepancy  has  been  offered  but  a 
tentative  explanation  Is  as  follows:  in  Long's  experiments  the  body  was 
moved  under  gravity  in  a  finite  cylinder  of  fluid  and  it  is  possible  that  the 
incipient  motion  of  the  fluid  in  ,  in  line  with  the  theory,  caused  a 
cavity  at  a  external  boundary  of  the  fluid,  the  filling  up  of  which  greatly 
distuibed  the  fluid  properties  in  .  It  would  be  of  Interest  to  examine 
the  flow  if  the  body  is  fixed  and  a  rotating  fluid  forced  past  it  to  see  whether 
the  fluid  still  closes  up  in  .  From  the  point  of  view  of  the  present  paper 
Long's  experlme  it,  while  making  us  optimistic  that  the  wake  in  C_  would 
be  obtained  experimentally,  must  make  us  cautious  about  accepting  the 
solution  with  =  -1 .  It  seems  in  fact  that  further  theoretical  progress 
depends  on  obtaining  experimental  evidence  of  the  flow  properties  at  large 


values  of  R 
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PART  U 

APPROXIMATELY  AUGNED  MAGNETIC  FIELD 


5.  Formulation 

Hera  we  suppose  that  the  motion  is  steady  but  that  the  imposed 
magnetic  field  is  inclined  at  a  small  angle  a  to  the  x-axls.  Further 
the  body  is  thin,  symne  trie  and  symmetrically  disposed  to  the  main 
stream  so  that  in  ( 3.  2) 


f^(x)  =  -f_(x)  . 

We  assume  that  the  perturijatlons  in  the  velocity  and  magnetic  fields 
are  small  and  write  for  their  components  relative  to  the  x,  y,  z  axes 


=  V^(1  + 


y’ 


0), 


ft  =  H^{cos  a  + 


sin  a  + 


h  , 
y’ 


0) 


(5.1) 


and  on  substituting  into  (2.1)  (2.  2)  a  set  of  linear  equations  is  obtained 
after  neglecting  squares  and  products  of  small  quantities.  It  may  then  be 
shown  that  can  be  expressed  in  terms  of  two  functions  i)j  viz: 


V  =  +  p[cos  a  +  sin  a  ] 

X  g^2  8x8y  gy2J 

2  2  2 

V  =  I'-l - P  [cos  a  +  sin  a  ] 

y  ax8y  ‘  2  8x8y  ^ 

OX 


1 

. 

j 


(5.2) 
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^2,  ^2,  \ 

h  B  COS  a  — ^  +  sin  or  ^  , 

X  g^2  8x8y  8x8y  ’ 


h  rcosa^  -  slnai^  . 

y  ®^®y  ax^  8x" 


(5.3) 


where  ^  Is  harmonic  and 


-  pRfcos^or^  +  2  cosofsina—^  +sin^a^  V  (5.4) 
^dx  8y  /  8x  \  8x  y  8y  / 


It  follows  that 


P“P» 


PV, 


go 


«  .0  + 

ax^ 


P(-cos 


.-di. 

8x8y 


+  sin  a  ) 
dx 


(5.5) 


The  equation  (5.4)  satisfied  by  ^  is  superficially  somewhat  different 
from  ( 2. 11)  but  they  are  reconcilable  if  regarded  as  particular  cases  of  a  more 
general  equation  in  which  a  ^  0  and  v|>  depends  on  t .  We  shall  see  however 
that  the  non-uniqueness  obtained  in  Part  I  is  absent  in  the  present  investigation. 
The  reason  for  this  difference  is  not  entirely  clear  and  further  work  on  this  aspect 
of  the  problem  would  seem  desirable. 

As  in  earlier  work  (10)  we  make  use  of  Fourier  transforms  to  get  a 
general  solution  of  t|j  in  a  form  convenient  to  our  purposes.  Write 


#360 


-29- 


♦(X,  y)=-^  /  (5.6) 

•  00 

00 

J  (5.7) 

•  00 

in  the  half  plane  y>0  where  '1'^  are  arbitrary  functions.  From  the 
differential  equations  satisfied  by  ^ ,  )|i 

=  +  (5.8) 

2  2  2  2  2  3 

[iu  PR  sin  a]x^  "  2ipRw  cos  a  sin  a  “  R  -  “  PR  cos  a  -  iw  =0  (5.  9) 

and  to  ensure  convergence  of  the  integrals  | ,  x^.  >nust  have  positive  real 
parts  on  the  real  axis  of  w  .  So  far  as  i  is  concerned  this  means  that  the 
upper  half  plane  re  u  >  0  must  be  cut  along  the  positive  imaginary  axis, 
the  lower  half  plane  re  w  <  0  must  be  cut  along  the  negative  imaginary 
axis  and  |  chosen  to  be  real  and  positive  on  the  real  axis  of  w .  The 
specification  of  x^  requires  a  little  more  care.  From  (5.  9) 
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J7 


when  «  is  large  x.^  =* -  and  it  follows  that  “Jto’  must  be  «|w| 

T  W 

being  positive  and  negative  with  u .  Hence  the  u  plane  must  be  cut  at 
infinity  both  above  and  below  the  real  axis  and  the  positive  sign  chosen 
in  (5. 10)  .  For  convenience  the  cuts  may  be  made  along  the  imaginary 
axis  and  therefore  they  extend  Inwards  towards  the  origin  as  far  as  the 
first  zeroes  met  of  the  expression  under  the  square  root  sign  in  (5. 10) 

1.  e.  as  far  as  Is^,  -Is^  where  s^,  -s^  are  the  positive  and 


negative  roots  of 


-  R(p-l)s  -  sln^  a  =  0  . 


From  these  points  Inwards  to  the  origin  the  cuts  are  absent  until  the  next 
zeroes  are  reached  and  since  these  both  occur  at  the  origin  the  factor 
47  may  be  replaced  by  u .  Thus 


itolpRcos  a  sin  a  +  -i<»>R(P-l)  ■¥  PR^sin^  g) 


iu  +  PR  sin  a 


If  R  =  « 


(5.11) 


(5.12) 


which  meai.s  that  the  corresponding  form  for  is 
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F^(x-c^y),  y>0  (5.13) 

where  Is  an  arbitrary  function,  In  agreement  with  earlier  work  (5). 

This  solution  represents  an  Alfv6n  wave  extending  to  Infinity  In  a  direction 
between  the  magnetic  field  and  the  downstream  directions.  If  a  s  0 ,  1.  e. 
the  magnetic  field  Is  aligned 

-  1«R(P-1)  (5.14) 

the  cut  extending  over  the  whole  Imaginary  axis  except  for  the  strip 

0  <  l|u|  <  1rIp-I|  .  Comparison  of  (5.13)  with  (5.15)  Indicates  that 

the  case  a  :  0  Is  not  a  straight  forward  limit  of  the  solution  when  a  #  0  . 

The  singular  nature  of  the  limit  Is  also  brought  out  by  studylr.g  the  double 

2  2 

limit  procedure  a-*0,  R  «  holding  Ror  =0  fixed.  The  expression 
for  Is  then 


^  =  -  -—j  {  pe  +  -  lu)(p-l) } 

R*  i«  +  pe 


(5.15) 


showing  how  sensitive  the  form  for  Is  to  the  value  of  0 . 

Parallel  arguments  to  the  above  may  be  used  to  show  that  when 


y  <  0 
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j  “  iwx  +  I  y 

«t>(x,  y)  =2^  J  •_(“)«  du 

•  00 

► 


.  «  i«x  +  X  y 

4'(x,  y)  =-^  J  ♦_(“)e  ■  dw  J 


(5. 16) 


where  ft  ,  ft  are  arbitrary  functions  and 


X  = 


iw 


!<*>  +  PR  sin  a 


-PR  cos  a  sin  a 


f  2 

+  Vw  - 


iu)R(p  -1)  +  pR^sin^a 


(5. 17) 


If  R  =  oo 

X  =  “7“  (  -  cos  or  +  p"^ )  =  iwc 
sin  a 


(5.18) 


leading  to  a  corresponding  form  for  <\i , 

F_(x  +  c_y)  (5.19) 

which  is  different  from  (5.13).  If  a  =  0  however 

X_  =  [w^-  lwR(P -1)  ]*  (5.20) 

which  is  Identical  with  ( 5. 14)  .  Thus  ifp>l,  aV^O,  R  =  oo  the  solutions 
for  k)j  is  y  >  0 ,  y  <  0  correspond  to  wakes  roughly  pointing  downstream 
and  upstream .  But  if  P>1,  a  =  o,  R<<»  the  two  solutions  are  of  similar 
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2  2 

type.  Finally  lfa-*0,  Rar  =0 

-4  =  ■  (-pe-t-N/pe^-  Mp-i)  )  (5.21) 

R^  lu  +  pe‘ 

which  is  different  from  (5.15)  for  all  6  >  0,  but  not  if  0  =  0. 

In  consequence  we  can  expect  the  wake  structure  of  the  flow  when 
or  =  0  to  be  of  a  different  kind  from  the  wake  structure  when  or  ^  0  ,  even 

-4 

if  O’  R  ^  so  that  in  some  sense  the  limit  a  -»  0  is  singular.  In  the  next 

2 

section  the  properties  of  the  solution  when  a  «  1 ,  R  »  1  but  aR  ~  1 
will  be  explored:  it  will  be  found  that  the  singular  nature  of  the  limit  is 
confined  to  the  wake  structure,  the  form  on  the  body  when  a  =  0  being 
deducible  from  the  form  when  a  4  0. 
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6.  Approximately  Aligned  Flald;  Symmetric  Body 

In  view  of  the  difference  between  the  wakes  for  p  >  1  predicted 

in  [4]  when  the  magnetic  field  is  oblique  and  ir  =  oo,  and  Lary's  theory, 

it  is  of  interest  to  examine  the  properties  of  the  solution  when  R  is  large 

and  a ,  the  angle  between  the  undisturbed  velocity  and  magnetic  fields,  is 

small.  It  will  be  supposed,  following  the  discussion  in  §5  that  the  orders 

2  2 

of  magnitude  are  such  that  Ra  =  0  is  of  order  unity.  Let 


[|^]=P(x),  [||]=Q(x),  [|^]=S(x),  [|i]=T(x)  (6.1) 

where  the  symbol  [  ]  denotes  the  leap  in  the  function  as  y  increases 
from  0-  to  0+  .  If  we  denote  the  transforms  of  P,  Q,  S,  T  by  bars. 


,  -_L.  9.  .  -  P  Q 

+  ■  21«  ■  2$  ’  -  2| 


(6.  2) 


^  _ T  ^  _ I_ 

*+"MX^  +  X.)  X^  +  X.  ’  -~’i"(X_+X^)  "x^+x 


(6.3) 


these  formulae  holding  for  general  values  of  a,  R.  We  now  specialize  to 

small  values  of  a  and  then,  from  (5.3),  since  h  ,  h  are  continuous 

X  y 


at  y  =  0 
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-  —  -  ^ 
P  +  aQ  +  T  =0 

- 

a  P  -  Q  +  S  =0 

4 


Further  If  we  write 


Vy  =  F{w)  sgn  y  +  Fj(w) 


(6.4) 


(6.5) 


at  y  =  0  where  F(u)  is  the  Fourier  transform  of  a  function  equal  to  f|(x) 
if  |x|<  1  and  is  otherwise  zero,  while  F^(u)  is  the  Fourier  transform  of  a 
function  which  is  arbitrary  in  |xi  <  1  and  zero  in  |x|<  1  ,  (5.2)  yields 


F(u)  =  ili-Sd-p) 


(6.6) 


rj(«) 


-•|eT  +  s 


«Px.x^ 

x^+x_ 


i‘^(X_  -  X^) 
2(X^.  +X_) 


(6.7) 


using  (6.4)  and  retaining  only  leading  terms.  From  (6.6)  we  have 


S'(x)  =-r^  f|(x)  |xl<l  (6.8) 

and  otherwise  is  zero,  in  agreement  with  the  corresponding  result  in  (  5) . 

Finally  since  the  pressure  is  continuous  on  that  part  of  the  x-axis  which 
does  not  include  the  body  we  have 
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iw(  P  -  apS  +  PT)  =  2G(u))  ( 6.  9) 

where  G(w)  Is  the  transform  of  a  function  which  is  zero  outside  |x|  =  1  . 

There  is  now  available  sufficient  information  to  determine  all  unknown 
functions  In  terms  of  f^(x)  apart  from  an  Indetermlnancy  associated  with 
circulation  and  usually  removed  by  a  suitable  Kutta  condition.  If  in  particular 
R»l,  ar  «  1  and  0  is  of  order  unity  x^»  X  are  given  by  (5.15), 

(5.21)  and  since  from  (6.4)(6.9)  T'(x)  is  non -zero  only  if  |x|  <  1  , 

(6.  7)  may  be  interpreted  as  an  integral  equation 


T'(x»)dx« 
X  -  x' 


-2pe(Up) 


/ 


X 


S’(x')dx' 

— 

(x'-x)* 


exp  -  { 


(x'-x)} 


(6. 10) 


giving  T'(x)  in  terms  of  S’(x)  apart  from  a  complementary  function  of  the 
2-4 

form  C(l-x  )  ‘  where  C  is  an  arbitrary  constant  presumable  determined 
by  a  suitable  Kutta  condition.  It  is  noted  that  from  (6.10)  T'  and  S'  are 
of  the  same  order  of  magnitude  and  that  T*  0  as  6  0  in  agreement  with 

Lary  (5). 

We  are  now  in  a  position  to  compute  the  force  on  the  body  and  to 
comment  on  the  wake  structure.  Restricting  attention  to  p  >  1  the  drag 


on  the  body  is 
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+1 

f 

af|(x)dx ( A +  A p  )dx 


(6.11) 


where  Ap^  are  the  excess  pressures  on  the  top  (y  =  0-f)  and  bottom 
surfaces  of  the  body.  Using  the  formula  (5.  5)  for  the  pressure  in  terms 
of  and  (6.2),  (6.3)  it  is  fairly  easy  to  show  that  the  Fourier 

transform  of  A  p^  +  A  p  is 


2PpV^R*F(u,) 

(pe^-Mp.l))2 


(6. 12) 


whence 


2PpV«9  lf;(x')dx'  2 

Ap  +  Ap_  = - 1  J  - T-  exp  -  (x'  -X)} 

a{w(p-l)}*  X  (x'-x)* 


(6.13) 


and  ( 6. 1 1)  may  now  be  evaluated  in  any  specific  case.  At  large  values  of 
6  (6.13)  reduces  to  2f^(x)/a  in  agreement  with  (4)  and  as  6 -»  0  , 

(6.13)  reduces  to 


(6. 14) 
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in  agreement  with  Lary  ( 5) .  From  the  formula  for  the  excess  pressure  it 
is  possible  to  deduce  a  condition  for  the  validity  of  the  assumption  that 
the  linearized  equation  and  boundary  conditions  may  be  used,  namely 

Max{R^,  ^}«j-  (6.15) 

where  c  is  a  measure  of  the  thickness  of  the  body. 

Let  us  now  look  at  the  structure  of  the  wake  when  P  >  1  ,  focusing 
attention  on  the  downstream  wake,  present  if  6  =  <»  but  absent  if  6  =  0. 

The  upstream  wake  and  the  case  P  <  1  can  be  examined  by  similar  methods. 

The  downstream  wake  is  given  by  8^i|//8x8y  for  it  makes  the  dominant 
contribution  to  p,  v^  in  the  range  of  values  of  R,  or  of  interest  here.  The 
wake  only  appears  if  y  >  0  when 


4  go 

8'>1)  _  .  f 

8xdy  “  "  2ir 


F(<o)do* 

{pe^-  1«(P-1)}^ 


exp [ lux 


Y  {pe  +  Vpe^-  iw(p-l)  }  .  (6. 16) 
w-ipe 


As  6  go  this  expression  reduces  to 


- 1  J  F(w)du  exp  iw  {x  -  c  y} 

ZttoP  *  -  00 


=  ,  f;(x.cy) 

ap2  ^  ^ 


(6. 17) 
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If  |x  -  c  y|  <  1  and  Is  otherwise  zero,  where  c  is  defined  in  (5.12). 

This  result  is  in  agreement  with  (4)  and  shows  that  when  6  is  large  the 

wake  is  controlled  by  the  exponential  term.  At  a  finite  value  of  6 ,  x  >  1  , 

and  taking  t^(x)  to  be  a  polynomial,  the  contour  of  integration  in  (6. 16) 

2 

may  be  deformed  into  a  finite  closed  curve  C  surrounding  u  =  ipo  . 

As  6  -•  0  this  contour  can  become  a  point  circle  at  the  origin  so  that 
the  integral  vanishes.  When  0  <<  1  only  the  behaviour  of  the  integrand 
near  w  =  0  is  relevant  and  since  f^(l)  =  “U  the  limit  of  F{u)/iuj  as 
cj-»  0  exists  and  we  write  it  as  A.  Consequently  (6.16)  is  equivalent  to 
writing 


^  Y) 


if  0  <<  1  where  X  =  0  x,  Y  =  0R*y  and 


(6. 18) 


^  (X,  Y)=-~  /  - — - 1  exp  {isX  (P  +  N/P  -  is(P-l)  }  . 

C  {p-ls(p-l)}*  ®  ^ 

If  Y  is  small,  ^  may  be  expanded  as  a  Maclaurln  series  in  Y  v«4iose 
coefficients  are  functions  of  X;  the  leading  terms  are 

Zp^e"^^  +  2p\e‘^^(PX- 3)  +  O(Y^)  . 


(6.19) 
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If  Y  is  large  ^  may  be  evaluated  by  the  method  of  steepest  descent 
and  we  have 

«cp .  {iSzziil ) 

T  (2irY®)*  2t  Y 


where  t  s  1  4  .  Thus  the  strength  of  the  wake  falls  off  rapidly  as  6  **  0 

from  ( 6. 18)  and  in  addition  becomes  much  more  diffuse  X,  Y  remaining 

small  for  a  wide  range  of  values  of  x,  y .  Nevertheless  the  characteristic 

wake  structure, albeit  in  a  very  weakened  form,  remains  for  all  6  >  0  .  This 

follows  from  ( 6.  20)  which  shows  that  the  wake  is  centred  on  the  line 

I 

X  s  tY  i.  e.  X  =  c^y  and  spreads  out  like  Y*  far  downstream  from  the 
body.  It  is  noted  that  ( 6.  20)  represents  the  wake  structure  when  Y  is 
large  for  all  non- zero  0  including  6  =  <»  and  consequently  shows  how 
the  structure  depends  on  0 .  For  example  when  6  is  large  the  thickness 
of  the  wake  ~  y^/6^  as  would  be  expected  since  the  wake  has  a  constant 
finite  thickness  when  6  -  «o  . 
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